Fixing Einstein's equations 



IFP-UNC-528 
kr-qc/990102l 



Arlen Anderson and James W. York, Jr. 
Dept. Physics and Astronomy, Univ. North Carolina, Chapel Hill NC 27599-3255 

(Jan. 6, 1999) 



Einstein's theory of general relativity has not only proven to be physically accurate [|T|, 
Q\ ', it also sets a standard for mathematical beauty and elegance — geometrically. When viewed 
as a dynamical system of equations for evolving initial data, however, these equations have 
a serious flaw: they cannot be proven to be well-posed (except in special coordinates 



That is, they do not produce unique solutions that depend smoothly on the initial data. To 
remedy this failing, there has been widespread interest recently in reformulating Einstein's 
theory as a hyperbolic system of differential equations p|-pO|] . The physical and geometrical 
content of the original theory remain unchanged, but dynamical evolution is made sound. 
Here we present a new hyperbolic formulation that is strikingly close to the space-plus-time 
( "3+1" ) form of Einstein's original equations. Indeed, the familiarity of its constituents make 
the existence of this formulation all the more unexpected. This is the most economical first- 
order symmetrizable hyperbolic formulation presently known to us that has only physical 
characteristic speeds, either zero or the speed of light, for all (non-matter) variables. It also 
serves as a foundation for unifying previous proposals. 

The source of the imperfection in Einstein's theory lies in the fact that it is a constrained 
theory. Physical initial data cannot be freely specified, yet even infinitesimally perturbed 
'. data that violate the physical constraints can lead to results so wildly divergent that they 
spoil the desired smooth dependence on initial data. This is particularly troublesome in 
numerical evolution where such violations are unavoidable. The lack of well-posedness is also 
^ ■ a serious problem when addressing such a basic question as the global nonlinear stability of 
flat Minkowski spacetime. In fact, the proof of such stability |H| employs the hyperbolic 
wave equation of |7j which we discuss below. A well-posed formulation of Einstein's equations 
would also seem to be an essential starting point for the conventional approach to quantum 
gravity in which one first quantizes the (unconstrained) classical theory and then imposes 
the constraints. 

The desire to simulate numerically the full nonlinear evolution of Einstein's equations 
in three dimensions, such as in the collision of two black holes (cf., e.g., H22|), has moti- 



vated much of the recent effort on hyperbolic formulations: well-posed underlying equations 
make stable numerical evolution much more likely than otherwise would be the case, and 
formulations cast in first-order symmetrizable form are especially suited to numerical imple- 
mentation. In addition, physical characteristic speeds make it easy to impose good boundary 
conditions, crucial to a successful numerical scheme. To amplify this point, Einstein's equa- 
tions contain many unphysical ("gauge") variables among its unknowns, and a priori they 
can travel at any speed. A formulation with only physical characteristic speeds has signifi- 
cant advantages because no explicit separation of physical and unphysical degrees of freedom 
is required. The physical and unphysical variables propagate at the same speeds and there- 
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fore satisfy boundary conditions on the same characteristic surfaces. This is particularly 
important, for example, at the horizon of a black hole, which is a characteristic boundary 
for physical variables but not for unphysical ones, unless the latter propagate at the speed 
of light. 

The following system of thirty equations will be shown to be symmetrizable hyperbolic 

= d o g ij + 2NK ij , (1) 

Rij = -N-%Kij + R% - N- 1 Vi VjN(a, g) + KK i3 - 2K ik K k 3 , (2) 

2g ij R k0 = d T kij + d 3 (NK kl ) + d t (NK kj ) - d k (NK l3 ) - 2 9ij NV m (K km - g km K) (3) 

(notation elaborated below). This form suggests the name "Einstein-Christoffel system." It 
is convenient for pedagogical reasons to replace the third equation by the equivalent equation 

*9k(iRj)o = doGkij + d k {2NK i3 ) - Ag k(l NV m (K j)m - g j)m K) (4) 

The indices i, j, k, m run over the spatial indices 1, 2, 3, and the Einstein convention of 
summing over repeated indices is used. A short-hand notation to indicate symmetrization, 
A^Bj) = (l/2)(A i Bj + AjBi), is used in the previous equation. 

To establish notation, assume that spacetime has topology £ x R with metric given in 
the foliation-adapted basis, 

ds 2 = -N 2 (dt) 2 + g^dx* + (3 i dt)(dx j + (3 j dt). (5) 

Here, N(a,g) is the lapse scalar, and (3 l (x,t) is the spatial shift vector, freely specifiable 
on the spacelike slices t = constant. The lapse iV is determined through iV = ag 1 ^ 2 , where 
a(x, i) is a freely specified "slicing" density (of weight —1) and g = det gi 3 is the determinant 
of the spatial metric g i3 . The spatial derivatives of the metric are denoted by 

Qkij = d k g i:j . (6) 

This is a subtle element, as it will transpire that while this relation is imposed initially, it 
may not hold for the evolved quantities, as discussed below. The spatial Christoffel symbols 
in this metric are, with T m ij = g mk T ki j, 

Fkij(G) = (l/2)(Q jki + Q ikj - Q kij ). (7) 

To focus attention on Q kij , the Christoffel symbols will not be used here as independent 
variables, though they could be, but only as a compact notation for this expression in terms 
of Qkij] the explicit functional dependence "((/)" is intended to reflect this choice. Finally, 
Kij denotes the extrinsic curvature of the slice S, and K = K k k is its trace. 

The derivative is the spatial covariant derivative operator in S. The derivative 
do = d t — £p, where d t = d/dt and £ p is the Lie derivative along the shift vector (3 in a 
t = constant slice, is the natural time derivative for evolving time-dependent spatial tensors. 
It is the extension to tensors of the (non-coordinate) basis vector do = d t —(3 k d k (d k = d/dx k ) 
that is normal to the slice S. Note that while [<9 , dj] = d dj — djd = (djf3 k )d k ^ 0, 
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[d ,d k ]=0. 



(8) 



On the left-hand side of (Q) and (|3|) or (|), and Rj are spacetime Ricci curvature 
tensors and are to be replaced by their appropriate expressions in terms of matter stress 
tensors from Einstein's equations, G^ v = R^ v — {\/2)g^ v R = 87rT^. Here, fi, v run from 
to 3, R = R^n, and T^ v is the matter stress tensor. Rfj is the spatial Ricci curvature 
tensor of the spacelike slice S. It is essential to manipulate the standard form of Rij = 
dkT k ij—djT k ik+T k m kT m ij—T k m jT m ik into a distinct but (initially) equivalent form, indicated 
below, and the superscript "(e)" reflects this change. 

Einsteinian initial data for the system ([!]), @, (|) are gr^, K^, Q ki j, specified on an initial 
slice S , and presumed to satisfy the Einstein constraints, G° = 8ttT° and R°k = 8iiT k . 



This system of initial constraints is well understood as a semilinear elliptic system p4| , |25 
A mathematically well-posed form of the twice-contracted Bianchi identities pE|JZT| shows 
that these initial-value constraints remain satisfied if the equations of motion are equivalent 
to Rij = 87r(T ij - (l/2) 9ij T%). 

That the system (JJ), (0), (f|) is hyperbolic is not obvious, but its content is easy to grasp. 
The first equation ([I]) is simply a definition of the extrinsic curvature K^. The second 
equation (0) is the 3+1 decomposition of the space-space components of the spacetime Ricci 
tensor. As such, these are the basic geometric ingredients of Einstein's original equations 
and of all 3+1 formulations of general relativity |2B|-^|. The remarkable fact is that the 
equation (f|) completes the first two into a symmetrizable hyperbolic system. The content 
of equation (|J) is also readily understood. 

If one applies d to @ and uses (§) and ([[]), one obtains the identity 

doOkij = -d k (2NKij). (9) 

This is the right hand side of @ aside from the 3+1 decomposition of Ag^R^Q. Ordinarily 
in Einstein's theory, 

Roj = -NV m {K rm -g jm K) (10) 

is a constraint — the "momentum" constraint — because it involves no time derivatives. What 
is special about (^) is that it makes the momentum constraint dynamical by combining it 
with the identity (§) involving a time-derivative. This defines a modified evolution of Qkij 
when the constraint is not satisfied, that is, when ([IT]) does not hold after R j is replaced 
by its matter expression. 

The identity ([D is closely related to metric compatibility of the connection. In a general 
spatial frame, a connection is metric compatible if and only if + = d^g^. Taking 
the time derivative of this condition and applying (||) gives 

2d T mk = -d k {2NKij) + Ag k{i NV rn (K ])m - g j)m K) + Ag k{i R j)0 . (11) 

Thus, if the momentum constraint were violated, metric compatibility of F would be lost in 
the evolution of (0). While 2T^ k = Q ki j always holds, the evolved Q k ij would no longer be 
the spatial derivative d k of the evolved g^. 

To motivate the system ([J), (|), ([|) further, consider two of its predecessors, the Einstein- 
Ricci formulation ||7|,|ll2l- |i4] , |i8|j and the Frittelli-Reula formulation [fT0|JT~6|] . The third-order 
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Einstein-Ricci system consists of ([[]) and a wave equation built from (Q) and ( JTOj) through 
the combination 



- Vi^-o - V^o = NnKij + Jij + S i:j . (12) 

It is called third-order because of the effective number of derivatives of gy in (|12|). Here, 
□ = —N~ 1 doN~ 1 do + V fe Vfc. is a nonlinear function of Kij, N, their first derivatives, 
and the second derivatives of N. Sij is a potentially troublesome term involving a second 
spatial derivative of K and a third derivative of N. The behavior of SV,- is tamed by using 
N = a(x,t)g 1 / 2 |7|,[T2"|j (or by imposing generalized harmonic slicing |18[ with a gauge source 



ID). Note that the use of a permits any time-slicing to be employed. 

This system can be put in first-order form by introducing new variables to represent 
the temporal and spatial derivatives of K%j and of N. Together with (fj) and the equations 
obtained by applying do to T k ij, and using a to eliminate N, a system of 66 equations is 
found. This system is manifestly spatially covariant, is expressed in 3+1 geometric variables, 
and has only physical characteristic speeds, either zero or the speed of light. 

One may wonder about the large number of equations and about a deeper meaning behind 
the combination in (JT2|). Regarding the number of equations, the Einstein-Ricci system is 
equivalent (for Einsteinian initial data) to the Einstein-Bianchi system |jn^,[2(J which also 
has 66 equations. There, it is evident that this number of equations is precisely that needed 
to incorporate the full Bianchi identities and to propagate the Riemann curvature tensor 
explicitly in a system having only physical characteristics (otherwise, cf. ||15||). 

The Frittelli-Reula system UlUyiB] , in contrast, has 30 equations, is expressed in non- 
covariant variables, and admits superluminal characteristic speeds for some (unphysical) 
degrees of freedom. It was constructed by proposing an energy norm built from the ex- 
trinsic curvature and spatial derivatives of the metric. Frittelli and Reula parametrize their 
construction and find a large family of hyperbolic systems with different characteristics, none 
wholly physical. Friedrich has observed that an equation for the metric constructed from 
(0) and (fj), while not of known hyperbolic type, has only physical characteristics A 



natural question is whether there are further thirty- variable hyperbolic systems and whether 
any has only physical characteristics. The Einstein-Christoffel system presented here is such 
a system, and it is easy from it to see how to extend the Frittelli-Reula construction. 

The Einstein-Ricci system has only physical characteristics, but its equations number 
more than twice those of the Frittelli-Reula system. A natural question is whether the 
third-order form can be put in first-order form to achieve a thirty variable system. The 
answer is yes. To see why this might be possible, consider the simple wave equation 

d\u - d 2 x u = 0, (13) 

This can be put in first order form in two ways. The easiest is to introduce the derivatives 
of the dependent variable as new variables. Introduce U = dtu and V = d x u to reach the 
system 

d t u = U, (14) 
d t U - d x V = 0, 
d t V - d x U = 0. 
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The last equation is an integrability condition reflecting the commutativity of the partial 
derivatives. This parallels the way that the first-order form of the Einstein-Ricci system was 
obtained from (|T2|). 

The second way to get to first order form is to pull apart the wave equation to obtain 
first order pieces 

d t u - d x v = 0, (15) 
dtv — d x u = 0. 

The first method is essentially one derivative higher. It necessarily involves more variables. 
Note that the wave equation ( |TB| ) is reconstructed from this system ( |T5| ) by taking a time 
derivative of the first equation and adding a spatial derivative of the second. This parallels 
the structure in flT2| ) that leads to the third-order Einstein-Ricci system. This encourages 
the speculation that a "pulled- apart" system analogous to (|Tj|) is possible. The obstacle 
however is that the momentum constraint as usually construed is not a dynamical equation, 
so the obvious "pulled-apart" system is not a hyperbolic system. The key idea is that adding 
a suitably chosen dynamical identity to the momentum constraint overcomes this obstacle 
and leads to a symmetrizable hyperbolic system. 

To begin, we work with Guj rather than F k ij- Focus on the derivatives of the Christoffel 
symbols contained in R^j — N'^-ViVjN. These are the essential terms from the standpoint 
of hyperbolicity. [Recall that N = ag 1 ^ 2 , so VjN = g 1 l 2 dja + f k j k (Q)g 1 ^ 2 a.} These terms 
are reorganized as follows: 

duT^^Q) - djT^iG) - d^ k jk (G) = (16) 

= 2 (®k[g km {Gjmi + Gimj — Gmij)) — dj(g rs Girs) — di(g rs Gj rS/ 

= \ (-dk{g km g m ij) + di\9 r '{Gr.3 ~ Gjrs)] + d^&rsi - Sirs)] 



+9 r g Sm {QirsQkmj + QjrsQkmi) g g Gkrs^Gjmi Gimj ) J ■ 



where d k (g Gjmi) = dj(g km Gkmi) + 9 ' g sm (GjrsGkmi - GkrsGjmi) has been used. Introducing 

fkij = l^Gkij ~ 9ki9 rS (Grsj — Gjrs) ~ 9kj9 TS (Grsi ~ Girs)} (17) 
~^9ki9 (J^rsj -Tjrs) -|- ~9kj9 (J^rsi -^irs) 

puts the leading derivatives of (0) in the familiar form 

Rii = -N^doK.j - d k f kij + Lo-y, (18) 

where l.o.y stands for lower order terms containing no derivatives of unknowns. They are 

l.o.ij = HKij - 2K ik K k j - a- l Vid 3 a - f k jk {G)ar l Via (19) 
+2T k mk (Q)T m ij(Q) — T k m j(Q)T m ik (g) — g kr g sm Gkrsfmij 

+^9 kr 9 sm (G jrs G kmi GirsGkmj) ~ 9 r 9 Sm Gkrs ( Gjmi Gimj ) ] • 
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Turn to consider (f|). From (|17|) , one computes 

9kiRjo + 9kjRio = -dofkij ~ dk^NKij) + l.o. kij . (20) 

The lower order terms are 

l.o. kij = NK ki g rs (Q rsj - Q jrs ) + NK kj g rs (Q rsi - Q irs ) (21) 
+g kl [K mj d m N - HdjN + Ng rs T k 3r {G)K sk - N(Q rsj - 2Q irs )K rs ] 
+g kj [K mt d m N - HdiN + Ng rs f k vr (G)K sk - N(G rm - 2G irs )K rs }. 

The subsystem (|T8|) and ([20| ) [completed by (P] is obviously symmetrizable hyperbolic 
because it has the familiar structure of a wave equation in first order form. [It is not 
rigorously symmetric hyperbolic because a metric is present to raise the spatial derivative 
index in ([18]).] It is also clear that to build a wave equation in from ([18]) and (f20D, 
one forms the combination in (|T2"D. This reveals the meaning behind this combination. The 
characteristic speed in the subsystem ([TJ), Q2l]| ) is the speed of light, so the extrinsic curvature 
and the connection propagate at the speed of light. From (||), the metric propagates at speed 
zero. 

It should be emphasized that the wave equation obtained from ([18]), (|20D is not exactly the 
same as the third-order Einstein- Ricci system. They differ by lower order terms proportional 
to constraints. This means that they will agree for Einsteinian initial data, but may disagree 
when the constraints are violated. 

The energy norm for the system ([!]), (|T8|), (PD|) is the integral over S of K^K^ + f kl ^f k ij, 
where / = g km g ir g^ s fmrs- When the energy norm is expressed in terms of Pij = Kij—gijK, 
P = P k k , and M ki j = —{\/2){G ki j — gijg rs G krs ), the result can be compared to the ansatzoi 



Frittelli-Reula ||16|| . Additional terms are present beyond those that they considered. Their 
energy norm ansatz can be generalized by adding arbitrary positive multiples of quadratic 
terms formed from traces of P^ and of M^ k . Similarly, their parametrization of Mij k can 
be extended using combinations of traces of Gijk, such as occur in f ki j. This produces a 
larger many-parameter family of symmetrizable hyperbolic systems equivalent to Einstein's 
equations. Some of these other systems also have only physical characteristics, for example, 
a multiple of the Hamiltonian constraint G°q can be added to (|2]) if a compensating multiple 
of the momentum constraint i? fc is added to (|3|) and (f|). 

The Einstein-Christoffel system discussed here is a well-posed system of 30 equations 
that has only physical characteristic speeds and can be expressed in 3+1 geometric variables. 
Spatial covariance of the system is not explicit, but present nonetheless. This formulation 
clarifies the relationships among Einstein's original equations and the Einstein-Ricci and 
Frittelli-Reula hyperbolic formulations. One can see further links to the Einstein-Bianchi, 



Friedrich P|6|,|8|,|T5|1 and Bona-Masso [§,0 formulations. It will be very interesting to im- 
plement this system numerically to see how it behaves. Having only physical characteristics 
should prove very useful when imposing boundary conditions on the horizon of a black hole 
since no information, physical or otherwise, can leave the black hole, and all information 
will enter the black hole in a physical way. 

Acknowledgement. AA thanks Mark Hannam for discussions. This work was supported 
in part by the NSF. 
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